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ABSTRACT 


This  report  studies  the  antenna  synthesis  characteristics  of  linear  phased  ar¬ 
rays  for  both  the  broadside  and  endfire  cases.  The  underlying  technique  is 
based  on  a  modification  of  the  conventional  Dolph-Chebyshev  side  lobe  ta¬ 
pering  technique.  Some  of  the  parameters  that  are  paramount  in  the  design 
of  these  arrays  are  analysed,  such  as  the  number  of  elements  N,  the  phase 
factor,  side  lobe  levels,  directivity,  radiation  pattern,  element  separation  and 
impact  upon  SNR.  The  results  are  compared  and  the  optimised  configuration 
is  investigated  for  each  case. 
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Phased  Array  Analysis  Using  a  Modified  Chebyshev 

Approach 


EXECUTIVE  SUMMARY 


The  work  presented  in  this  report  is  partially  aimed  towards  the  Microwave  Radar  Branch’s 
Microwave  RAdar  Test  Environment  (MIRATE)  project  which  has  the  objectives  of  al¬ 
lowing  microwave  radar  researchers  and  engineers  to  build  skills  and  first  hand  knowledge 
of  phased  array  radars  and  to  generally  foster  R  &  D  that  ultimately  will  be  used  for  the 
analysis  of  more  complex  phased  array  radar  systems  such  as  the  one  anticipated  to  go  into 
the  new  Air  Warfare  Destroyer  (Project  SEA  4000).  Specifically,  one  component  under 
MIRATE  is  the  construction  of  an  operational  experimental  phased  array  radar  (XPAR), 
which  as  a  first  step,  comprises  the  development  of  a  16-element  linear  array.  This  report 
presents  theoretical  results  that  can  be  tested  experimentally  by  the  XPAR  program  and  it 
is  hoped  that  other  methods  and  techniques  will  be  studied  and  implemented  throughout 
various  stages  of  XPAR. 

The  approach  here  is  to  consider  an  arbitrary  number  of  linear  arrays  and  study  the  var¬ 
ious  design  specifications  that  might  surface  for  such  a  system  using  a  method  based  on 
conventional  Dolph-Chebyshev  tapering  techniques.  By  using  a  modified  Chebyshev  ap¬ 
proach  (MC)  we  investigate  amongst  other  things  the  side  lobe  level,  directivity,  radiation 
pattern,  element  spacing,  relation  to  the  number  of  elements  and  so  on.  The  MC  method 
has  superior  convergence  properties  for  large  arrays  given  that  the  window  coefficients 
of  the  radiation  pattern  can  be  determined  via  the  use  of  matrix  algebra.  At  the  same 
time  this  avoids  the  use  of  complicated  polynomial  equations  that  are  computationally 
intense.  Results  presented  are  for  both  broadside  and  endfire  arrays  and  comparison 
shows  interesting  behaviour  depending  on  the  selection  of  either  an  even  or  odd  number 
of  elements. 
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1  Introduction 


In  this  report  we  will  consider  an  analysis  that  involves  a  variation  of  the  technique  used 
in  conventional  Chebyshev  arrays.  A  key  aspect  of  Chebyshev  arrays  is  that  they  produce 
equal  side  lobe  levels  (SLL)  for  a  given  radiation  pattern.  However,  Chebyshev  arrays 
suffer  from  directivity  saturation  when  the  number  of  radiating  elements  becomes  large 
[1].  Thus  there  is  a  need  to  find  a  method  to  analyse  large  Chebyshev  arrays  in  an  efficient 
manner.  Furthermore,  the  synthesis  of  the  radiation  pattern  is  difficult  because  for  every 
radiation  pattern  a  new  complicated  polynomial  series  needs  to  be  found.  It  is  of  no 
surprise  that  researchers  have  been  looking  for  alternatives  in  order  to  overcome  these 
inefficiencies  [2].  What  we  shall  call  the  Modified  Chebyshev  (MC)  method,  is  a  technique 
that  is  based  on  the  synthesis  of  arrays  using  the  zeros  of  conventional  Chebyshev  arrays 
repeatedly,  see  for  example  [3].  More  specifically,  the  idea  behind  the  MC  formulation  is 
that  it  makes  no  direct  use  of  Chebyshev  polynomials.  The  array  factor  is  calculated  in 
terms  of  cosine  and  hyperbolic  cosine  functions  respectively,  while  a  system  of  equations 
for  the  excitation  amplitudes  is  obtained  from  the  zeros  of  the  array  factor.  The  solution 
of  the  system  of  equations  for  the  excitation  amplitudes  is  done  by  allowing  one  of  those 
excitation  amplitudes  to  be  the  independent  variable,  as  we  shall  see  later.  The  procedure 
allows  for  the  solution  of  any  number  of  radiating  elements  N,  for  both  an  even  or  odd 
configuration  of  arrays  respectively.  Lee  [4]  has  used  the  zeros  of  the  array  factor  before  in 
order  to  study  a  different  formulation  of  Chebyshev  arrays,  however  the  method  still  uses 
a  polynomial  representation  for  the  array  factor,  which  does  not  generalise  and  becomes 
mathematically  involved  when  we  consider  a  large  number  of  elements  or  large  arrays. 
The  problem  of  beam  synthesis  for  example,  involving  an  unknown  number  of  elements 
but  specified  SLL  and  beamwidth,  is  tackled  traditionally  by  determining  the  minimum 
number  of  elements  that  are  required  to  meet  the  design  parameters  through  a  trial  and 
error  process.  In  what  follows,  we  will  examine  a  process  that  allows  the  minimum  number 
of  elements  to  be  obtained  in  a  single  step  without  resorting  to  an  iterative  process,  hence 
giving  a  more  flexible  and  efficient  design  strategy.  Furthermore  we  assume  that  among 
the  SLL,  half-power  beamwidth,  and  number  of  elements,  two  quantities  are  specified  and 
the  third  quantity  as  well  as  the  excitation  amplitudes  ^  are  sought.  In  addition,  we  will 
also  consider  both  broadside  and  endfire  arrays  as  well  as  the  optimum  spacing  between 
the  radiating  elements. 


2  The  Array  Factor 


We  consider  a  linear  array  of  N  equally  spaced  elements  with  uniform  excitations  and 
symmetrical  distributions  for  the  magnitude  amplitudes  in  the  array,  as  well  as  a  constant 
inter-element  phase  shift  a  of  radiating  electromagnetic  waves  at  angle  9  from  the  plane 
of  the  array-refer  to  Figure  1.  For  an  even  and  odd  number  of  elements  we  can  express 
the  array  factor  as. 


^These  are  obtained  relative  to  the  amplitude  of  one  element  that  is  taken  as  the  independent  variable 
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Figure  1:  Linear  array  with  elements  radiating  at  angle  6  and  phase-shift  a.  The  inter¬ 
element  separation  is  d. 


n  r  /  s  1 

/(V')  =  2  E  «mCos  ■,N  =  2n, 


ffif)  =  ao  +  2  E  «m  cos  [mf)]  ;  iV  =  2n  +  1,  (2) 

m=l 

where  if  =  pdcos{6)  +  a,  /3  =  27r/A;  A  is  the  wavelength,  d  is  the  element  spacing,  a  is 
the  phase  factor,  0  is  the  angle  measured  from  the  line  of  the  array,  Om  is  the  magnitude 
of  the  amplitude  for  the  element  on  either  side  of  the  array  midpoint  and  oq  denotes 
the  amplitude  of  the  centre  element  when  N  is  odd.  Given  the  fact  that  Chebyshev 
polynomials  satisfy  the  relationships  r„(cosx)  =  cos(nx)  for  |T„|  <  1  and  T„cos(/ix)  = 
cosh(nx)  for  >  1,  the  array  factor  in  (1)  or  (2)  is  written  as 


\fW\  =  \fW  = 


C  cos  (^- 
C cosh  ^ 


(N-l), 

9  ‘ 

'(N-l) 

2 


\m\<c 

\fm>c, 


where  C  is  a  positive  constant  coefficient  and  is  related  to  by 


cos  (IJ 
cosh  (I) 


\fm<c 

\fm>c. 


In  (4)  we  note  that  7  is  a  real  coefficient  that  can  be  determined  when  the  design  specifica¬ 
tions  are  considered.  By  close  examination  of  (3),  we  find  that  /{ip)  assumes  its  maximum 
value  when  ip  =  0.  Moreover,  this  corresponds  to  ip  =  ipo  =  2cosh“^(7).  In  most  cases 
of  interest,  we  can  assume  a  knowledge  of  the  number  of  radiating  elements  N  and  the 
maximum-to-side  lobe  level  ratio  R  ^  .  Since  we  assume  that  the  side  lobes  have  the  same 
magnitude  (see  (3)),  we  can  write  down  the  ratio  R  as 


(iV-1), 


R  =  cosh 


By  transposing  (5)  and  solving  for  7  we  obtain: 


ipo  )  =  cosh  (iV  —  1)  cosh 


7  =  cosh  loge(-R  +  VR'^  -1) 


^The  side  lobe  level  here  is:  —20logioiR)  in  dB. 
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Figure  2:  (i)  For  an  even  number  of  elements  N,  the  amplitude  distributions  are  obtained 
by  normalising  all  other  amplitudes  by  the  edge  element  an-  The  ealeulations  in  Tables 
1-2  refer  to  this  eonfiguration  for  even  N,  where  the  symmetry  point  is  taken  to  be  at  a 
distanee  of  d/2  from  the  element  with  amplitude  ai-  (ii)  the  same  as  before  exeept  we  now 
eonsider  an  odd  number  of  elements  N  with  the  symmetry  oeeuring  at  the  element  with 
amplitude  qq  at  a  distanee  of  d  from  ai . 


We  are  now  in  a  position  to  calculate  7  when  we  know  N  and  the  half-power  beamwidth 
HPBW  or  alternatively  we  can  do  the  same  if  we  know  R  and  the  HPBW  as  we  shall  see 
later.  Finally,  if  we  require  the  array  factor  to  be  normalised  to  unity,  we  can  determine 
the  coefficient  C  =  1/R. 


2.1  The  Excitation  Amplitudes 


A  system  of  equations  that  govern  the  magnitude  of  the  radiating  elements  can  be  con¬ 
structed  from  the  zeros  of  the  array  factor  in  terms  of  matrices.  These  zeros  are  so¬ 
lutions  of  cos[(Ai  —  1)^/2]  =  0,  and  when  rearranging  for  we  obtain  the  solutions 
ip  =  [(4m  ±  l)7r]/(A^  —  1),  with  m  being  an  integer.  The  fact  that  the  array  factor  is 
symmetrical  about  ■0  =  tt  allows  us  to  consider  only  the  zeros  in  the  range  D  <  <  vr-see 

Figure  2  for  the  distribution  of  the  elements.  The  zeros  are  now  obtained  from. 


_  (2m  —  l)7r 
“  iV-  1 


(7) 


for  m  =  1,2,3,  ...,M,  while  M  =  {N  —  2)/2  for  N  even  and  M  =  {N  —  l)/2  for  N  odd 
number  of  radiating  elements  respectively.  More  explicitly  we  obtain  the  zeros  in  terms  of 
Tp  by 


'Ipm  =  2  COS 


'1 

—  COS 

L7 


/  (2m  —  l)7r  Y 

V  2(iV-l)  )\  ’ 


(8) 
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Table  1:  Amplitude  magnitudes  for  equally  spaeed  Modified  Chebyshev  arrays.  Note  that 
the  amplitudes  of  the  edge  elements  are  always  unity  and  all  other  amplitudes  are  nor¬ 
malised  by  the  former,  see  Figure  2.  The  Table  below  shows  values  sueh  that  the  first  row 
eorresponds  to  ai  (even)  or  qq  (odd),  the  seeond  row  eorresponds  to  02,....,  a^-i,  0^  =  1. 


N 

\SLL\  =  -10  dB 

\SLL\  =  -20  dB 

\SLL\  =  -30  dB 

\SLL\  =  -40  dB 

3 

1.0390 

1.6364 

1.8774 

1.9604 

1 

1 

1 

1 

4 

0.8794 

1.7357 

2.3309 

2.6688 

1 

1 

1 

1 

5 

0.7975 

1.9319 

3.1397 

4.1448 

0.7248 

1.6085 

2.4123 

3.0131 

1 

1 

1 

1 

6 

0.6808 

1.8499 

3.3828 

4.9891 

0.6071 

1.4369 

2.3129 

3.0853 

1 

1 

1 

1 

7 

0.6102 

1.8387 

3.7846 

6.2731 

0.5864 

1.6837 

2.3071 

5.2678 

0.5191 

1.2764 

2.1507 

3.0071 

1 

1 

1 

1 

8 

0.5413 

1.7244 

3.8136 

6.8448 

0.5103 

1.5091 

3.0965 

5.1982 

0.4519 

1.1386 

1.9783 

2.8605 

1 

1 

1 

1 

9 

0.4923 

1.6627 

3.9565 

7.6989 

0.4813 

1.5800 

3.6516 

6.9168 

0.4494 

1.3503 

2.8462 

4.9516 

0.3995 

1.0231 

1.8158 

2.6901 

1 

1 

1 

1 

10 

0.4463 

1.5585 

3.8830 

7.9837 

0.4306 

1.4360 

3.4095 

6.6982 

0.4003 

1.2125 

2.5986 

4.6319 

0.3576 

0.9264 

1.6695 

2.5182 

1 

1 

1 

1 

11 

0.4113 

1.4907 

3.8985 

8.4813 

0.4054 

1.4421 

3.6982 

7.8954 

0.3880 

1.3036 

3.1457 

6.3341 

0.3601 

1.0949 

2.3702 

4.2952 

0.3235 

0.8450 

1.5400 

2.3546 

1 

1 

1 

1 

12 

0.3787 

1.4031 

3.7865 

8.5669 

0.3697 

1.3277 

3.4657 

7.5913 

0.3521 

1.1860 

2.8885 

5.9124 

0.3269 

0.9948 

2.1659 

3.9692 

0.2952 

0.7759 

1.4262 

2.2035 

1 

1 

1 

1 
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Table  2:  Continued  from  Table  1.  Amplitude  magnitudes  for  equally  spaeed  Modified 
Chebyshev  arrays.  Note  that  the  amplitudes  of  the  edge  elements  are  always  unity  and  all 
other  amplitudes  are  normalised  by  the  former,  see  Figure  2.  The  Table  below  shows  values 
sueh  that  the  first  row  eorresponds  to  ai  (even)  or  oq  (odd),  the  seeond  row  eorresponds 
to  02, - ,  On—\,  On  —  !• 


N  \SLL\  =  -10  dB 

\SLL\  =  -20  dB 

\SLL\  =  -30  dB 

\SLL\  =  -40  dB 

13  0.3528 

1.3408 

3.7438 

8.8242 

0.3492 

1.3101 

3.6071 

8.3860 

0.3387 

1.2218 

3.2208 

7.1776 

0.3217 

1.0830 

2.6498 

5.4821 

0.2990 

0.9093 

1.9856 

3.6660 

0.2713 

0.7167 

1.3261 

2.0658 

1 

1 

1 

1 

14  0.3286 

1.2683 

3.6224 

8.7888 

0.3229 

1.2189 

3.3979 

8.0481 

0.3119 

1.1245 

2.9818 

6.7204 

0.2958 

0.9931 

2.4335 

5.0683 

0.2753 

0.8358 

1.8273 

3.3899 

0.2510 

0.6655 

1.2379 

1.9409 

1 

1 

1 

1 

15  0.3086 

1.2138 

3.5549 

8.8944 

0.3064 

1.1932 

3.4582 

8.5628 

0.2996 

1.1330 

3.1807 

7.6287 

0.2885 

1.0378 

2.7574 

6.2585 

0.2735 

0.9146 

2.2401 

4.6829 

0.2549 

0.7724 

1.6883 

3.1411 

0.2335 

0.6209 

1.1598 

1.8281 

1 

1 

1 

1 

where  m  =  1,2,3,  We  make  the  observation  that  the  V'm’s  are  zeros  of  /(V^)  in  (1) 

or  (2),  that  is. 


/(V'm)  =  0,  m  =  1,2,3,. ..,M. 


(9) 


In  fact  (9)  gives  a  system  of  M  equations  with  M  +  1  unknowns.  These  unknowns  to 
be  solved  are  the  amplitudes  of  the  radiating  elements  om,  one  of  which  is  chosen  as  the 
independent  variable.  In  the  case  of  an  even  number  of  elements,  the  independent  variable 
is  chosen  to  be  om+i-  We  can  form  matrices  that  allow  us  to  solve  for  any  number  of 
array  elements  N .  From  (1)  we  have  the  form 


n 

f{'4)k)  =  2  X!  «mCOS 

m=l 


(10) 
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for  k  =  1,2,3,  The  solution  of  (10)  for  the  amplitudes  Um  can  be  calculated  from 

the  form 


r  n 

T 

r  ^  ^  ^  1 

T' 

ai  a2  •  •  •  aM 

=  — OM+l 

7i  72  •••  7m 

cos[(M  +  l/2)V'i]  cos[(M  +  1/2)V'2]  •  •  •  cos[(M  +  1/2)^^^] 


(11) 


which  can  be  abbreviated  to. 


a  =  — OM+i 


(12) 


where  T  is  the  transpose.  We  will  examine  (12)  and  therefore  define  the  matrices  A,  C  and 
the  vectors  7,  by  solving  (10)  via  (12)  for  /c  =  1,2.  Furthermore,  we  will  consider  N  =  6 
elements  for  the  array,  and  because  we  are  dealing  with  an  even  number  of  elements  we 
have  M  =  2.  From  this  we  know  that  the  independent  variable  we  need  is  given  by  om+i 
or  03.  Expanding  (10)  we  obtain 


fi'ipi)  =  2ai  cos 


2 


+  202  cos 


+  203  cos 


(f), 


and 


f{'ip2)  =  2ai  cos  +  2a2  cos 

Now  from  (9)  we  know  that  the  expansions  in  (13)  and  (14)  are  zero,  ie. 


(13) 


(14) 


(15) 


Equating  (13)  and  (14)  to  zero  and  transposing  the  independent  amplitude  and  its  coeffi¬ 
cient  to  one  side,  it  is  easy  to  show  that  we  obtain  the  matrices: 


cos 


cos 


( 


Ol 


=  -03 


\  02 


COS 


cos 


(16) 


We  recall  that  03  in  (16)  is  the  independent  variable  and  has  been  transposed  to  the  right 
side.  We  are  of  course  interested  in  solving  (16)  for  the  amplitudes  and  so  we  have 


( 


Ol 


cos 


=  -03 


\  02 


A 


cos 


cos 


cos 


cos 


(17) 


The  solutions  for  ai  and  02  satisfy  the  array  factor  (1)  or  (10),  for  an  even  number  of 
radiating  elements  N .  If  we  let 


-.T 
a  = 


C  = 


iT 


Ol  02 


COS 


.T 

A  = 


(+) 


cos 


(+) 

iT 


(18) 


7i  72 
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71  =  _  cos  -  cos  ^  _ 

72  =  [  cos  (^)  cos  (^)  ]  ’ 

we  recover  (12)  for  the  case  k  =  1,2.  The  procedure  we  have  studied  so  far  means  that 
we  can  generalise  the  results  above  to  any  number  of  radiating  elements,  and  it  is  easy  to 
extend  to  the  general  form; 


=  —aM+i[A  ]■ 


cos[(M  + 
cos[(M  +  ^)'lp2] 

cos[(M  +  ^)'iPm] 


where  the  matrix  A  has  the  general  form, 


A  = 


cos 


(*,) 


cos  (Af  + 
cos  (M  +  \)^ 


Here  M  =  {N  —  2)/2  for  an  even  number  of  array  elements  N.  In  a  similar  way  we  can 
now  consider  the  case  where  N  is  an  odd  number  of  elements.  For  example,  if  Ai  =  5,  we 
obtain  from  M  =  {N  —  l)/2  that  M  =  2  and  we  also  recall  that  the  independent  variable 
is  chosen  to  be  oq  •  The  array  factor  (2)  can  be  written  as 

n 

fi4>k)  =  ao  +  2  E  dm  cos[mV’fc],  (22) 

m=l 

so  if  we  seek  solutions  for  k  =  1,  and  2  for  instance  we  obtain  the  expansions 

/(V’l)  =  flo  +  2ai  cos  (V’l)  +  202  cos  {2ipi)  (23) 

and 

/('02)  =  ao  +  2ai  cos  (V’2)  +  202  cos  (2V’2)  •  (24) 

From  (15)  we  find  once  again  that  (23)  and  (24)  can  be  written  as 

/  cos(V^i)  cos(2V^i)  W  /  00  \ 

i  cos  (^^2)  cos  (2V^2)  \a2  \  «o  /  ’ 


which  finally  gives  the  solutions 


--2«o 


cos  (V’l)  cos  (2 V’l) 
cos  (V’2)  cos  (2 V’2) 


As  we  have  done  so  before,  we  can  now  generalise  for  odd  N  by  writing  the  generalised 
form  as 


-Vo  [s 


-1  1 
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or  in  the  vector-matrix  form 


a 


I, 


-.T 

where  the  matrix  B  becomes 


B 


/  cos(V’i) 

cos  (V’2) 


cos  (2^^!) 
cos  (2V'2) 


\  cos(V’m)  cos(2V’m) 


cos  (MV’i) 
cos  {M1P2) 


V 


cos(MV’m)  / 


(28) 


(29) 


The  solutions  for  am  satisfy  the  array  factor  as  given  by  (2)  for  an  odd  number  of  ele¬ 
ments  N.  Table  1  shows  the  calculated  amplitude  coefficients  for  even  and  odd  radiating 
elements  using  the  technique  discussed  above.  Also,  the  method  investigated  so  far  hints 
at  the  possibility  of  obtaining  the  amplitudes  via  a  series.  Indeed  this  is  the  case,  such  a 
series  solution  is  based  on  the  work  by  Barbiere  [5]  but  for  large  arrays  there  are  conver¬ 
gence  problems.  For  small  arrays  it  gives  excellent  results  that  allow  comparison  of  the 
amplitudes  for  an  even  or  odd  number  of  elements.  We  first  define  the  variable  zq  to  be: 


(i?  +  \/i?2  _  (^-1)  +  (i?  - 


(30) 


where  R  once  again  represents  the  ratio  of  the  side  lobe  level.  For  an  even  number  of 
elements,  M  =  N/2,  and  the  amplitudes  are  given  by 


Cln 


E 


q=n 


{2M  -  l){q  +  M  -  2)\  .M-q  2q-l 

{q  —  n)\{q  +  n  —  1)\{M  —  q)\ 


(31) 


Alternatively,  for  an  odd  number  of  elements  M  =  (N  —  l)/2,  and  thus  we  obtain  the 
amplitudes  by  the  series  expansion. 


Cln 


M-1 


q=n 


i{Q 


2M{q  +  M  -2)1 _ 

n)l{q  +  n  —  2)!(M  —  g  -|-  1)! 


,\M-q+l2(q-l) 


(32) 


with  the  additional  condition  that  is  given  by 


2,  n  =  1 
1,  n  /  1 


(33) 


such  that  n  =  1,  2, 3, ....  It  is  now  a  matter  of  expanding  the  even  and  odd  solutions  of  the 
array  factors  respectively, 

M 

AF^even)  =  COS  [(2n  -  1)«]  (34) 

n=l 


and 


8 


M+l 

^p{odd)  _  ^  anCOS  [(2n  —  l)u], 

n=l 


(35) 
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Array  Factor 


Figure  3:  Left  Plot;  The  array  factor  as  a  function  of  u  for  \SLL\  =  —10  dB.  The 
curve  in  red  represents  N  =  17  (odd)  number  of  elements  while  the  curve  in  blue  is  for 
N  =  20  (even).  Right  Plot;  Polar  plot  for  the  parameters  used  on  the  left. 


U 

Figure  4-  Left  Plot;  The  array  factor  as  a  function  of  u  for  |S'LL|  =  —10  dB.  The 
curve  in  red  represents  =  11  (odd)  number  of  elements  while  the  curve  in  blue  is  for 
N  =  10  (even).  Right  Plot;  Polar  plot  for  the  parameters  used  on  the  left. 


Array  Factor 


curve  in  red  represents  N  =  17  (odd)  number  of  elements  while  the  curve  in  blue  is  for 
N  =  20  (even).  Right  Plot;  Polar  plot  for  the  parameters  used  on  the  left. 
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Table  3:  Polynomial  approximation  for  the  HPBW  using  the  optimum  spaeing  between 
the  radiating  elements  N  for  a  broadside  array.  The  left  hand  eolumn  represents  the  side 


lobe  ratio,  while  the  right  side  eolumn  is  the  approximation  to  the  half-power  beamwidth 
degrees  for  any  N  =  3, ...,  20. 

\SLL\  (dB) 

H  PBW  {degrees) 

-10 

-2.97864  X  10“‘^iV'''  -h  2.65506  x  10“‘^iV^  -  9.89469  x  10“^iV^ 
-h0.1998571V^  -  2.3685iV3  +  16.59961V2  _  65.1958iV  -h  120.406 

-20 

-3.06897  X  10“*^ W''  -h  2.75572  x  lO-^A^*^  -  1.03612  x  10“^^^ 
-h0.211604iV^  -  2.54401iV3  -h  18.1839Af2  -  73.4929iV  -h  141.942 

-30 

-2.50017  X  10-‘^iV''  -h  2.26546  x  lO-^iV^  -  8.61964  x  10-^iV^ 
-h0.17886iV^  -  2.19797iV3  -h  16.2084Ar2  _  68.6181V  -h  142.699 

-40 

1.199797  X  10“^iV^  -  1.5909  x  10“^Af^  -h  5.18349  x  10“^iV^ 
-0.889623iV3  8.60381iV2   46.0883iV  -h  119.759 

where  we  define  u  =  ^cos(0),  and  d  as  being  the  separation  between  the  radiating 
elements.  The  results  are  exact  to  those  obtained  by  the  MC  method  and  are  thus  both 
displayed  in  Table  1.  Since  the  variable  u  G  [—1,1],  corresponding  to  0  =  180  and  0 
respectively,  we  can  obtain  the  array  factor.  Generally,  using  this  type  of  method,  it 
appears  that  an  odd  number  of  radiating  elements  give  a  broader  beamwidth,  but  greater 
’directivity’  or  gain  than  an  even  number  of  elements.  Interestingly,  odd  arrays  can  give 
’better’  synthesis  with  less  number  of  elements  but  the  correct  choice  will  also  depend  on 
the  requirements,  eg,  where  the  nulls  appear.  Figures  3-5  show  the  radiation  pattern  for 
different  N  and  SLL.  An  interesting  result  arises  when  we  consider  the  behaviour  of  large 
arrays  as  a  function  of  the  SLL.  Figure  6  shows  a  plot  of  the  ’average’  of  the  radiation 
amplitudes  <  a  >  as  a  function  of  the  array  elements  N,  for  a  given  SLL.  In  particular 
for  a  -10  dB  SLL^,  we  can  investigate  the  contibutions  of  the  average  of  the  amplitudes 
radiated  by  the  array  elements  by  the  use  of  linear  regression  techniques.  We  define  the 
average  <  a  >  to  be 

i 

where  a*  are  the  individual  amplitudes  from  each  of  the  array  elements  and  we  consider  a 
linear  fit  to  the  data  in  order  to  obtain  the  expression 

<  a  >=  1.00474  -  0.04062551V  (37) 

where  N  is  the  number  of  elements  for  both  the  even  and  odd  cases.  Equation  (37)  allows 
us  to  solve  for  <  a  >  if  the  number  of  radiating  elements  are  given  or  alternatively  if  <  a  > 
is  given,  we  can  obtain  N.  Clearly  the  linearity  of  (37)  is  a  very  inaccurate  approximation 
and  a  better  fit  to  the  data  gives  the  result. 


<a>=  1.40296  -  0.144412iV  -h  0.00688^^  -  0.00012231V3.  (38) 

Once  again  (38)  allows  us  to  solve  for  <  a  >  or  N,  however  when  solving  (38)  for  N,  the 
cubic  equation  can  give  a  number  of  solutions,  some  of  which  will  be  imaginary.  In  order 

®The  SLL  =  —lOdB  example  is  used  here  in  order  to  demonstrate  the  method.  The  same  procedure 
therefore  applies  for  any  SLL. 
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Table  Polynomial  approximation  for  the  HPBW  using  the  optimum  spaeing  between 
the  radiating  elements  N  for  an  endfire  array.  The  left  hand  eolumn  represents  the  side 
lobe  ratio,  while  the  right  side  eolumn  is  the  approximation  to  the  half-power  beamwidth  in 
degrees  for  any  N  =  3, 20. 


\SLL\  (dB) 

H  PBW  {degrees) 

-10 

-2.69065  X  10-^ N'  -h  2.41777  x  lO-^A^*^  -  9.10298  x  lO-^iYH 
0.186397iV^  -  2.25303Ar3  +  l6.3049Af2  _  68.3012iV  -h  160.484 

-20 

-2.3982  X  10“‘^A^'''  -h  2.16965  x  10“‘^iV*^  -  8.24329  x  10“3iV5 
-h0.170909iV4  -  2.10251A^3  15,6131V2  _  68.0538iV  -h  170.946 

-30 

-1.74512  X  10-‘^iV''  -h  1.59827  x  10-‘^iV^  -  6.16729  x  lO-^iV^ 
-h0.1304931V^  -  1.650611V3  -h  12.7594iV2  -  59.1631V  -h  165.512 

-40 

-1.07953  X  10“‘^A1'''  -h  1.00824  x  10““^ 77*^  -  3.98904  x  lO-^AT^ 

-h8. 71996  X  lO-^AT^   i.i5209iV3  9.4574iiv2  -  47.82771V  -h  153.68 

SLL=-10  dB 


Figure  6:  The  solid  eurve  is  obtained  from  (38),  while  the  points  represent  the  number 
of  elements  N  for  eaeh  average  power  <  a  >. 
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to  obtain  the  required  solution,  (37)  is  solved  to  obtain  a  value  Nq,  such  that  when  solving 
(38)  it  is  used  as  the  initial  guess  or  initial  condition  thereby  improving  considerably  the 
accuracy  of  (38).  We  can  analyse  other  values  of  the  side  lobe  ratio  in  a  similar  way. 
Furthermore,  the  use  of  regression  techniques  means  that  the  HPBW  (in  degrees)  for  both 
broadside  and  endfire  arrays  can  be  obtained  for  a  given  SLL  if  any  number  N  for  the 
elements  is  substituted.  This  is  shown  in  Table  3  and  Table  4  respectively. 


3  The  Power  Aperture  Product  of  a  Weighted 

Array 


From  previous  results  we  have  seen  that  generally  an  odd  number  of  elements,  while  giving 
a  slightly  broader  beamwidth,  has  greater  directivity  (or  gain)  than  an  even  number  of 
elements.  Interestingly,  when  we  consider  such  arrays  in  terms  of  aperture  power  and 
efficiency  relative  to  an  unweighted  array  the  opposite  of  what  we  expect  occurs.  For 
instance  as  we  will  discuss  below,  the  optimum  efficiency  of  the  array  is  greatest  for  an 
even  number  of  arrays  even  though  we  have  a  reduced  gain  as  a  result.  Therefore  the  type 
of  performance  we  require  will  determine  the  kind  of  trade-offs  we  are  prepared  to  accept 
in  antenna  synthesis.  The  power  aperture  product  of  a  weighted  array  is  given  by: 


Parr  ay  ^array 


■  N 

E 

Li=l 


r,^ 


Pelement  ^element 


(39) 


where  we  define  Parray  as  the  average  power  of  the  array,  Aarray  as  the  effective  area  of  the 
array,  Peiement  as  the  average  power  of  the  unweighted  element  (assumed  identical  across 
the  array),  G element  as  the  gain  of  the  unweighted  element  (assumed  identical  across  the 
array),  at^i  are  the  normalised  amplitude  weights  applied  on  transmission  to  element  i, 
Or-y  are  the  normalised  amplitude  weights  applied  on  reception  to  element  i.  Where  the 
element  normalisation  is  such  that  max{ai)  =  1,  it  reflects  the  practical  issue  that  the 
antenna  weighting  is  normally  achieved  by  attenuating  an  element.  For  an  unweighted 
antenna  dt,i  =  1  and  dr^i  =  1,  reducing  the  previous  result  to  the  well  known 


Parray  Aarray 


—  Peiement  Aelement 


(40) 


The  relative  efficiency  of  the  weighted  antenna  to  the  unweighted  antenna  is  then  given 
by: 


e  = 


1 


■  N 

^  )  ^t,idr^i 
.i=l 


PelementG element  1 


p. 


element  G  element 


■  N 

^  )  dt,idr^i 
.i=l 


(41) 


We  firstly  consider  the  performance  of  arrays  with  Chebyshev  weighting  applied  only  once, 
so  that  either  dt^i  =  1  or  dr^i  =  1.  Figures  7-9  show  the  power  aperture  product  and  the 
efficiency  of  the  single  weighted  array  for  differing  numbers  of  elements.  Each  of  the  Figures 
differ  by  the  maximum  number  of  elements  shown  in  order  to  highlight  the  bahaviour  of  the 
arrays  as  the  number  of  elements  increases.  What  is  apparent  from  the  Figures  is  that  for 
any  given  weighting  the  total  power  aperture  product  increases  as  the  number  of  elements 
is  increased,  but  that  beyond  a  certain  critical  number  of  elements  the  power  aperture  rate 
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Total  array  power  Rx  vs  number  of  elements  Efficiency  vs  number  of  elements 


Figure  7;  Left  Plot;  The  total  transmitted  aperture  power  is  shown  vs  a  relatively  small 
number  of  array  elements  N  at  different  sidelobe  levels.  Right  Plot;  The  effieieney  of 
the  aperture  power  of  arrays  is  shown  vs  the  number  of  elements. 


of  increase  is  sharply  reduced,  and  the  overall  power  aperture  efficiency  sharply  declines. 
The  reason  for  this  decrease  is  due  to  the  fact  that  the  Chebyshev  weights  must  transfer 
power  from  the  main  lobe  to  the  sidelobes  in  order  to  maintain  the  required  sidelobe  ratio. 

We  will  now  consider  if  for  a  given  two  way  side  lobe  level  it  is  more  efficient  to  weight  both 
transmission  and  reception  or  just  one  (usually  reception  because  of  interference  rejection 
considerations).  Figure  10  shows  the  total  power  aperture  product  for  two  cases.  In  the 
first  case,  20dB  Chebyshev  weighting  is  applied  to  both  transmission  and  reception,  and  in 
the  second  case  40dB  Chebyshev  weighting  is  applied  to  the  reception  signal  only.  Figure 
10  broadly  shows  that  for  small  arrays  it  is  better  to  apply  weights  on  both  transmission 
and  reception,  but  for  larger  arrays  it  is  more  efficient  to  apply  Chebyshev  weights  to 
reception  only^. 


4  The  Impact  on  the  Signal  to  Noise  Ratio  of  a 

Weighted  Array 


Whilst  the  signal  levels  are  impacted  by  the  weighting  function  applied,  so  are  the  noise 
levels  (at  least  on  receive),  and  so  we  must  consider  the  overall  impact  on  the  resultant 
signal  to  noise  ratio  (SNR). 

If  we  assume  that  the  noise  at  each  receiver  is  independent  and  of  the  same  mean  power 
level,  then  the  expected  value  of  the  total  noise  after  weighting  will  be: 


■^For  the  very  smallest  arrays  there  may  be  a  reversal  of  this  result,  as  shown  in  Fignre  10 
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Total  array  power  Rx  vs  number  of  elements  Efficiency  vs  number  of  elements 


Figure  8:  Left  Plot;  The  total  transmitted  aperture  power  is  shown  vs  N  =  100  array 
elements  at  different  sidelobe  levels.  Right  Plot;  The  effeieney  of  the  aperture  power  of 
the  arrays  is  shown  vs  the  number  of  elements. 


Efficiency  vs  number  of  elements 


Figure  9:  Left  Plot;  The  total  transmitted  aperture  power  is  shown  vs  a  relatively  large 
number  of  array  elements  at  different  sidelobe  levels.  Right  Plot;  The  effeieney  of  the 
aperture  power  of  the  arrays  is  shown  vs  the  number  of  elements. 
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Efficiency  vs  number  of  elements 


Total  array  power  Rx  vs  number  of  elements 


Figure  10:  Left  Plot;  The  total  transmitted  aperture  power  is  shown  vs  the  number  of 
elements  at  different  sidelobe  levels  for  a  two-way  transmit  and  reeeive  SLL.  Comparison 
is  made  of  the  power  for  an  unweighted  transmition  but  weighted  reeeption  of  SLL=-40  dB 
to  that  of  a  transmit  and  reeeive  mode  of  SLL=-20  dB.  Right  Plot;  The  effieieney  of  the 
aperture  power  of  the  arrays  is  shown  vs  the  number  of  elements. 


n^array  — 
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.i=l 


"^element 


(42) 


where  Uarray  IS  the  average  noise  power  of  the  array,  neiement  is  the  average  noise  power 
of  the  unweighted  element  (assumed  identical  across  the  array),  dr^i  are  the  normalised 
amplitude  weights  applied  on  reception  to  element  i.  The  signal  to  noise  ratio  is: 
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where 

{SNR),i,^,^t  =  (44) 

'^element 

An  array  that  is  unweighted  on  both  transmit  and  receive  gives: 

{S  N  R)un—^gigfited  —  N  {S  N  K)  element 
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SNR  ...  ,/SNR  Chebyshev  weights,  Rx  only 

weighted  ur-weighted  '  a  <  i 


Figure  11:  Impact  upon  SNR  of  Chebyshev  weights  applied  on  reception  only. 


Hence  the  impact  upon  the  SNR  of  weighting  the  antenna  is  given  by: 


{S  N  R)  weighted 

'  N 

i=l 

2 

{SN  R)  un— weighted 

N 

'  N 

E  al, 

4  =  1 

(46) 


or 
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y~) 


i=l 


N 


N 

E- 

i=l 


-  {S N weighted 


(47) 


The  impact  of  applying  Chebyschev  weights  only  on  receive  (that  is  at^i  =  1)  on  the  system 
SNR  is  shown  in  Figure  11,  with  the  loss  in  signal  to  noise  again  becoming  significant  once 
the  natural  side  lobe  levels  are  lower  than  the  Chebyschev  adjusted  level.  The  losses 
shown  here  broadly  match  those  reported  in  the  literature-see  for  example  Nathanson  [7] 
provided  the  array  is  sufficiently  small,  otherwise  the  literature  understates  the  losses  that 
are  experienced. 

Using  weights  on  both  transmission  and  reception  has  a  far  more  profound  impact  upon 
system  SNR,  as  shown  in  Figure  12,  with  a  significant  loss  in  SNR  resulting  from  the  loss 
in  power  aperture  efficiency  on  transmit  dominating  the  overall  SNR  performance.  This 
remains  the  case  where  the  overall  weighting  appears  nominally  equal  to  the  receive  only 
case  shown  in  Figure  13. 
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SNR  .  . .  .  /  SNR  .  u.  ^  -  Chebyshev  weights  applied  to  BOTH  Tx  &  Rx 

weighted  un-weighted  ‘  o  rr 


Figure  12:  Impact  upon  SNR  of  Chebyshev  weights  applied  both  on  transmit  and  reeep- 
tion. 


-  Chebyshev  weights,  Tx  8.  Rx  comparison 


Figure  13:  Impact  upon  SNR  of  Chebyshev  weights  applied  on  transmit  and  reception 
eompared  with  reception  only. 
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fx-Band  =  10  GHz 


Figure  I4:  The  optimised  spaeing  plotted  against  the  number  of  elements  for  an  X- 
Band  (X  =  0.03  m)  array.  The  eurves  starting  from  the  bottom  to  the  top  eorrespond  to 
|5LL|  =  —10,— 20,— 30  and  —40  dB’s  respeetively. 

5  Investigation  of  Beamwidth,  Side  Lobe  Level 
and  the  Number  of  Elements 


Another  important  parameter  of  interest  that  we  will  look  at  is  the  beamwidth.  We 
consider  9  =  Oh  such  that  it  corresponds  to  the  half-power  point  on  the  main  beam  of 
the  array  pattern  with  the  appropriate  phase  factor  a  chosen  according  to  whether  we 
have  broadside  or  endfire  arrays.  There  are  two  half-power  points  for  a  broadside  array 
(a  =  0),  one  exists  on  the  left  and  the  other  on  the  right  side  of  0  =  7r/2  respectively.  For 
the  case  Oh  <  71/2,  the  half-power  beamwidth  for  broadside  arrays  is  HPBW  =  vr  —  20h- 
On  the  other  hand,  for  endfire  arrays  (a  =  Pfdd),  there  exists  only  one  solution  for  Oh 
so  that  HPBW  =  Oh  when  a  =  —fdd  and  HPBW  =  11  —  0  when  a  =  +f5d.  Here  we 
will  consider  a  =  —fdd  for  the  endfire  case  without  loss  of  generality  [3] ,  and  note  that  for 
d  =  A/2,  Oh  =  HPBW/2.  We  define  if  =  ifn  at  0  =  Oh,  so  that  for  the  broadside  case 
'ipH  =  Pdcos{6H),  and  for  the  endfire  case  ifn  =  Pd[cos{0H)  —  !]•  From  (3)  we  have 


tpH 


2  cosh  ^ 


(48) 


We  recall  that  the  array  factor  assumes  its  maximum  value  at  =  V’O)  so  that  from 
equations  (3),  (5)  and  (6)  after  combining  we  obtain 


cos  (  1  cosh 


1 


cosh  ^  {R) 


N-1 


—  cosh 


=  0. 


(49) 


Equation  (49)  gives  the  relationship  between  the  half-power  beamwidth,  side  lobe  level 
and  the  number  of  elements.  More  specifically  if  two  of  these  are  given,  the  third  one  can 
be  calculated.  Thus,  for  the  case  where  N  and  HPBW  are  given,  we  can  determine  R. 
As  can  be  seen  from  (49)  however,  the  equation  involves  transcendental  functions  and  so 
the  solution  for  R  needs  to  be  found  numerically.  Any  number  of  solutions  can  satisfy 
(49)  so  that  in  order  to  pick  the  correct  solution  we  must  firstly  guess  a  ’root’  as  an  initial 


18 


DSTO-TR-1806 


fL-Band  =  1-2  GHz  fs-Band  =  3.0  GHz 


Figure  15:  The  optimised  spaeing  plotted  against  the  number  of  ele¬ 
ments  for  an  L-Band  (X  =  0.25  m)  array  on  the  left  and  an  S-Band 
fA  =  0.1  m)  array  on  the  right.  The  eurves  starting  from  the  bottom  to 
the  top  eorrespond  to  \SLL\  =  —10,  —20,  —30  and  —40  dB’s  respeetively. 


condition  before  proceeding.  This  means  that  (49)  needs  to  be  used  with  caution.  For 
example,  if  we  want  to  find  the  number  of  elements  N  for  a  given  R  and  HPBW,  we 
find  that  is  a  non-integer  number  which  physically  is  impractical.  The  correct  solution 
would  be 

iV  =  int  {A^o}  +  1,  (50) 

where  int  {A^o}  rneans  the  integer  part  of  the  solution  Nq  obtained  numerically  from  (49). 
To  further  highlight  the  method,  we  consider  an  endfire  array  with  \SLL\  =  —20  dB  and 
a  half-power  beamwidth  of  22.5®.  By  the  use  of  (49)  we  find  the  numerical  solution  to  be 
Nq  =  13.596  but  because  of  (50)  and  the  fact  that  we  need  to  have  HPBW  <  22.50®,  the 
actual  solution  is  Ai  =  14.  The  beamwidth  is  obtained  from 


HPBW  =  TT  —  2  cos  ^  )  for  broadside 

pd 


and 


HPBW  =  cos  ^  (  1  —  )  for  end  —  fire. 


pd 


where  from  (49)  we  transpose  to  get 

pH  =  2  cos” 


1  1 

cosh 

1  rosh“V  ^  1 

yv-l  ^osn  1^1 

1 

cosh 

[ Aril  cosh  ^(i?)] 

(51) 


(52) 


(53) 


For  endfire  arrays  where  the  inter-element  spacing  is  d  =  A/2,  we  multiply  the  right-hand 
side  of  (52)  by  the  factor  2.  Figures  14-15  show  results  for  frequency  bands  of  interest. 
For  selected  SLL,  the  optimum  configuration  is  obtained  by  the  given  optimum  element 
separation  and  number  of  elements.  In  fact  all  variables  can  be  calculated  simultaneously 
by  looking  at  the  graphs,  ie,  if  we  require  the  optimum  performance  for  the  array  then  for 
a  selected  optimum  element  separation  we  read  off  the  corresponding  number  of  elements 
and  SLL  and  so  on.  Similarly,  the  variation  of  the  HPBW  (in  degrees)  as  obtained  for 
optimum  element  separation  and  A/2  separation  for  both  broadside  and  endfire  arrays  is 
shown  in  Figures  16-20. 


19 


DSTO-TR-1806 


6  Optimising  the  Inter-Element  Separation 


We  are  interested  in  obtaining  the  distance  between  the  elements,  d,  such  that  for  a 
specified  SLL  and  given  number  of  elements  the  smallest  beamwidth  can  be  achieved. 
This  means  ’packing’  as  many  side  lobes  as  possible  within  the  radiation  pattern  of  the 
array.  In  fact  by  doing  so,  the  visible  range  in  the  V’  domain,  which  corresponds  to  the 
range  0  <  0  >  tt,  is  determined  in  such  a  way  as  to  include  all  side  lobes  on  both  sides  of 
the  main  beam  for  broadside  arrays  and  on  one  side  of  it  for  endfire  arrays.  A  portion  of 
the  grating  lobe  which  will  produce  a  side  lobe  with  a  magnitude  equal  to  the  specified 
SLL  is  also  included.  The  cutoff  point  on  the  grating  lobe  satisfies  the  relation  [3] : 


cosh 


[  2 

cos 

[  2  1 

(54) 


that  is  valid  only  if  ^  =  0.  The  solution  of  if)  which  corresponds  to  the  grating  lobe 
lies  in  the  range  tt  <  'ip  <  2tt.  In  this  range  however,  cos('0/2)  <  0  and  (4)  should  be 
modified  as  — 7cos('0/2)  =  cos(V;/2).  If  we  put  ?/)  =  0  in  the  previous  expression  we  obtain 
^p  =  2tt  —  2cos“^(1/7).  For  broadside  arrays,  xp  should  correspond  to  0  =  0,  that  is  when 
xp  =  fddopt,  where  dopt  is  the  optimum  spacing  between  the  elements.  As  a  result  we  have 


dopt  —  ^ 


cos  ^(-) 

2 _ lj_ 

TT 


broadside. 


(55) 


In  the  case  of  endfire  arrays  with  the  main  beam  in  the  6  =  tt  direction,  xp  should  corre¬ 
spond  to  0  =  0,  in  other  words  xp  =  2f3dopt,  resulting  in 


cos  H^)' 

TT 


end  —  fire. 


(56) 


Equation  (56)  is  also  valid  for  endfire  arrays  with  the  main  beam  in  the  0  =  0,  (a  =  —(3d), 
direction.  Optimisation  of  the  element  spacing  d  for  endfire  arrays  has  also  been  studied 
for  the  case  a  /  ±(3d,  see  [6] . 


7  Directivity 


In  this  section  we  will  consider  the  directivity  of  arrays  with  arbitrary  spacing  between 
the  radiating  elements.  Stegen  [7]  has  investigated  the  case  when  the  elements  have  a 
spacing  that  is  an  integer  multiple  of  half  a  wavelength.  Approximations  have  also  been 
obtained  for  the  directivity  of  a  large  number  of  elements  [8].  The  more  interesting  case 
of  an  analytical  representation  for  the  directivity  in  the  general  case  of  arbitrary  element 
spacing  and  arbitrary  number  of  elements  is  very  useful  [9].  In  what  follows,  we  will  obtain 
the  results  that  allow  us  to  calculate  the  directivity  of  arrays  as  a  function  of  the  number 
of  elements,  element  spacing,  and  side  lobe  level.  The  results  are  particularly  useful  when 
the  number  of  elements  are  not  large  and  the  element  spacing  is  not  an  integer  multiple 
of  half  a  wavelength.  We  begin  by  recalling  the  form  of  the  array  factor  for  a  Chebyshev 
array  which  consists  of  N  equally  spaced  isotropic  elements. 
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Broadside  Array  End-fire  Array 


Figure  16:  Comparison  of  the  half-power  beamwidth  (HPBW)  for  a 
broadside  array  (left)  and  an  endfire  array  (right)  with  optimum  sparing 
between  the  elements.  The  dashed  eurve’s  side  lobe  ratio  is  -10  dB,  while 
the  solid  eurve  is  that  for  a  -fO  dB  side  lobe  ratio. 


fW  =  ^  cos 


(A^  — l)cos  ^(tcosI-^)) 


when  I /('(/')  I  <  1/-R,  and 


/(V’)  =  ^  cosh 


(A^  — l)cosh  ^(7cos(-^)) 


(57) 


(58) 


when  |/('!/')|  >  1/R.  As  before,  if  =  (ddcosO  +  a,  d  is  the  element  spacing,  a  is  the  inter¬ 
element  phase  shift,  /3  =  27r/A,  A  is  the  wavelength  and  0  is  the  angle  measured  from  the 


line  that  is  perpend 
ratio  and  7  =  cosh 


icular  to  the  array.  Here  R  is  the  beam  maximum  to  side  lobe  level 
ln(i7  -|-  \/ R?  —  T)/{N  —  1)1 .  The  expression  for  the  directivity  can  be 


obtained  from  many  publications,  for  example  see  [10],  and  is  generally  denoted  as: 


D  = 
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TT  27r 

/  f  p{0)sm(e)ded<j> 
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f 
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(59) 


In  order  to  evaluate  the  integral  in  the  denominator  of  (59),  we  notice  that  the  expres¬ 
sion  /('0)  as  given  by  (57)  or  (58)  can  be  expanded  as  a  series  in  cos(nx)  or  cosh(nx) 
respectively.  By  collecting  terms  we  arrive  at  the  result  for  /^(V’): 
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where  M  =  N  —  1  and  xo  =  27cos('0/2).  By  substituting  (60)  into  (59)  and  using  the 
identity. 


/  = 2-="  ( ) 


2n  \  sin  [(2n  —  2k)x] 


(2n  -  2k) 


we  obtain  the  directivity  D  as 


D  = 


4R^ 


(61) 


(62) 
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N  N 

Figure  17:  The  half-power  beamwidth  (HPBW)  in  degrees  plotted  against  the  number 
of  elements  for  a  broadside  array  (left)  and  an  endfire  array  (right).  The  dashed  eurve 
represents  an  array  sparing  of  \ (2  while  the  solid  eurve  is  the  optimised  sparing.  The 
magnitude  of  the  side  lobe  level  is  -10  dB. 


Broadside  Array  End-fire  Array 


0  5  10  15  20  0  5  10  15  20 


N  N 

Figure  18:  The  half-power  beamwidth  (HPBW)  in  degrees  plotted  against  the  number 
of  elements  for  a  broadside  array  (left)  and  an  endfire  array  (right).  The  dashed  eurve 
represents  an  array  sparing  0/  A/2  while  the  solid  eurve  is  the  optimised  sparing.  The 
magnitude  of  the  side  lobe  level  is  -20  dB. 


Broadside  Array  End-fire  Array 


Figure  19:  The  half-power  beamwidth  (HPBW)  in  degrees  plotted  against  the  number 
of  elements  for  a  broadside  array  (left)  and  an  endfire  array  (right).  The  dashed  eurve 
represents  an  array  sparing  0/  A/2  while  the  solid  eurve  is  the  optimised  sparing.  The 
magnitude  of  the  side  lobe  level  is  -30  dB. 
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Broadside  Array  End-fire  Array 


Figure  20:  The  half-power  beamwidth  (HPBW )  in  degrees  plotted 
against  the  number  of  elements  for  a  broadside  array  (left)  and  an  end- 
fire  array  (right).  The  dashed  eurve  represents  an  array  spaeing  of  Xj^ 
while  the  solid  eurve  is  the  optimised  spaeing.  The  magnitude  of  the  side 
lobe  level  is  -fO  dB. 
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(64) 

(65) 


fid  (M  —  i  —  k) 

Equation  (64)  can  be  simplified  further  if  we  consider  a  broadside  array  (a  =  0)  or  endfire 
array  (a  =  Pfid)  so  that  for  both  of  these  cases  we  have 


C{i,k)  = 


sin  [v{M  —  i  —  k)fid] 
vfid{M  —  i  —  k) 


where 


V  =  1 
u  =  2 


broadside 
end  —  fire. 


(66) 


(67) 


Calculations  show  that  for  the  same  number  of  elements  and  side  lobe  level,  the  directivity 
of  a  broadside  array  with  arbitrary  element  spacing  d,  is  related  to  the  directivity  of  an 
endfire  array  with  an  element  spacing  of  d/2,  ie,  (see  Fig  21) 


k^broadside(.d)  Z^end— /ire(^^/2)  • 


(68) 


Furthermore,  when  d  is  a  multiple  integer  of  A/2  and  fid  =  mir,  where  m  is  an  integer,  we 
find  that  C,{i,k)  =  0.  For  this  latter  case,  r  in  (62)  becomes: 
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Thus  we  find  that  when  d  =  A/2,  the  directivity  of  broadside  and  endfire  arrays  is 
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Figure  21:  The  relation  between  the  SLL,  the  ratio  of  the  optimum 
spaeing  to  the  wavelength  (dopt/^)  and  the  number  of  elements  N.  The 
plot  on  the  left  is  for  a  broadside  array  while  the  plot  on  the  right  is 
for  an  endfire  array. 


X  X 

Figure  22:  Plots  of  the  direetivity  D  against  the  ratio  d/X  for  a  broad¬ 
side  array  (the  endfire  ease  differs  by  a  faetor  of  1/2).  The  number  of 
elements  that  have  been  ehosen  are  N  =  3,7, 12, 16, 19  and  20,  where 
N  =  3  is  the  lowest  eurve  and  N  =  20  is  the  highest  in  eaeh  plot.  The 
side  lobe  level  is:  -10  dB  for  the  left  plot  and  -20  dB  for  the  right 
plot. 
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Figure  23:  Plots  of  the  directivity  D  against  the  ratio  d/X  for  a  broad¬ 
side  array  (the  endfire  case  differs  by  a  factor  of  1/2).  The  number  of 
elements  that  have  been  chosen  are  N  =  3,7, 12, 16, 19  and  20,  where 
N  =  3  is  the  lowest  curve  and  N  =  20  is  the  highest  in  each  plot.  The 
side  lobe  level  is:  -30  dB  for  the  left  plot  and  -40  dB  for  the  right 
plot. 


equivalent.  For  a  specified  number  of  elements  N,  side  lobe  level  and  element  spacing  d/X, 
we  can  apply  (62)-(65)  to  calculate  the  directivity.  On  the  other  hand  if  d  is  a  integer 
multiple  of  A/2,  we  can  use  the  expression  for  r  as  given  by  (69).  Figures  (22-23)  show 
plots  of  the  directivity  against  the  ratio  d/X  for  broadside  arrays  with  the  endfire  curves 
differing  only  by  a  factor  of  1/2.  From  the  graphs  we  can  see  that  the  curves  ’peak’  for  given 
values  of  the  ratio  d/X.  The  element  spacing  d  (in  d/X)  corresponding  to  the  maximum 
of  the  directivity  is  not  equal  to  the  optimum  spacing,  ie,  d  /  dopt-  In  other  words,  we 
find  that  minimum  beamwidth  and  maximum  directivity  do  not  occur  simultaneously. 
The  results  presented  here  are  very  useful  when  we  are  dealing  with  situations  where  the 
element  spacing  is  different  from  integral  multiples  of  half  a  wavelength  and  the  number 
of  elements  is  not  too  large. 


8  Conclusion 


We  have  studied  arrays  based  on  a  modification  of  conventional  Chebyshev  arrays.  Cheby- 
shev  arrays  have  the  property  that  the  side  lobe  levels  are  the  same  but  are  hindered  by  the 
fact  that  the  directivity  ’saturates’,  especially  as  the  number  of  elements  increases.  The 
MC  method  also  has  the  same  side  lobe  levels  but  the  directivity  is  vastly  improved  even 
when  using  the  same  parameters  as  for  the  conventional  Chebyshev  method.  The  higher 
directivity  in  the  MC  method  is  achieved  by  finding  the  zeros  of  conventional  Chebyshev 
arrays  repeatedly.  However  the  maximum  directivity  is  achieved  when  the  spacing  between 
the  elements  is  slightly  smaller  than  the  optimum  element  spacing  in  the  case  of  the  MC 
method.  The  3  dB  beamwidth  narrows  or  broadens  depending  on  the  number  of  elements 
chosen  with  pronounced  differences  if  even  or  odd  elements  are  chosen.  All  parameters  for 
such  arrays  have  been  studied  and  comparisons  have  also  been  made  between  broadside 
and  endfire  arrays. 
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